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ABSTRACT 

Our XMM-Newton spectrum of the giant, high-redshift (z = 1.88) radio galaxy 6C 0905+39 
shows that it contains one of the most powerful, high-redshift, Compton-thick quasars known. 
Its spectrum is very hard above 2 keV. The steep XMM-Newton spectrum below that energy is 
shown to be due to extended emission from the radio bridge using Chandra data. The nucleus 
of 6C 0905+39 has a column density of 3.5 



-1.4 
-0.4 



x 10 24 cm 2 and absorption-corrected X-ray 



luminosity of 1.7_q X 10 erg s in the 2 — 10 keV band. A lower redshift active galaxy 
in the same field, SDSS J090808.36+394313.6, may also be Compton-thick. 



1 INTRODUCTION 



The powerful, classical double FRII dFanaroff & Rilevl 19741) radio 
galaxy 6C 0905+39 is unusual as it is one of the highest redshift 
(z = 1. 88) giant radio galaxie s known, spanning lllarcsec on 
the sky dLaw-Green et alJl 19951) corresponding to a projected size 
of 945 kpc in the cosmology assumed in this paper. Radio galax- 
ies have long been associated with obscured, (optically defined) 
Type 2 quasars, because they are considered to be the same ob- 
jects as radio -lou d quasars seen thro ugh an obscuring torus (e.g. 
lBarthel|[l989l and lAntonuccil 119931) . X-ray observations are par- 
ticularly efficient at detecting obscured AGN as they provide the 
observer with several diagnostics, such as the photoelectric cut- 
off, the iron line equivalent width (EW) and the X-ray to optical 
or [O III] A500 7 ratio (e.g. iBassani et alj [l 999 or, for a review, 
IComastril 12004). Sources with an obscuring column greater than 
1.5 x 10 24 cm -2 , at which the Thomson scattering optical depth 
reaches unity, are known as Compton-thick AGN. It is non-trivial 
to c alculate the column den sity of such sources as there is little 
(see lWilman & F abian 1999) direct emission below lOkeV mean- 
ing that analysis relies on interpreting the reflected spectrum at 
lower X-ray energies. Suzaku is starting to provide a window into 
the X-ray regime above 10 keV enabling direct detection of nuclear 
emission in the b rightest Compton-thick sources such as NGC 5728 
dComastri et alj |2007). Such sources are important since the X-ray 
background gives strong evidence for a significant populatio n of 
Compton-thick AGN (e.g. lComastrfe004IWorslev et all200H) . 

In this paper, we present recent XMM-Newton observa- 
tions of 6C 0905+39 (which is located at RA 09h08ml6.9s, Dec 
+39d43m26s in J2000). Throughout this paper, all errors are quoted 
at lcr unless otherwise stated and the assumed cosmology is 
H = 71kms _1 Mpc" 1 , fi = 1 and A = 0.73. One arc- 
second represents 8.518 kpc on the plane of the sky at the red- 
shift (z = 1.883 ± 0.003) of 6C 0905+39 and the Galactic ab- 



sorption along the line-of- sight is taken to be 1.91 x 10 20 cm 
dDickev & Lockmanlll99(il) . 



2 DATA REDUCTION 

XMM-Newton observed 6C 0905+39 on 2006 October 30 for 

58.0 ks with EPIC-pn and 61.9 ks with each of the EPIC-MOS 
cameras. The data were reduced using the standard pipeline, run- 
ning the S AS tools EPCHAIN and EMCHAIN for the EPIC-pn and 
both EPIC-MOS cameras respectively. The resulting files were fil- 
tered to remove flares, leaving 42.2 ks of good-time for the PN, and 

53.1 ks for both MOS 1 and MOS 2 when dead-time intervals are 
also taken into account. 

Spectra for the background (an area of sky free from sources 
near 6C 0905+39 on the same chip), the central source (within a 
15 arcsec circle) a nd the extended X -ray emission which is dis- 
cussed elsewhere dErlund et all 120071) were extracted separately 
for each instrument and stacked and fitted using XSPECv 11.3.2. 
Given the low number of counts in the MOS spectra, only the PN 
spectra were used for analysis of the core. 



3 RESULTS 

The X-ray emission extending along the radio axis over hundreds 
of kpc to the east an d west, which was fi rst detected by Chan- 
dra and presented in iBlundell et af] (2006), is clearly visible in 
the XMM-Newton data (see Fig. [TJ. In total there are 2566 counts 
(1783 background) in the source, 2125 (1624 background) in the 
extended emission and 441 counts (159 background) in the core 
in the 0.2 - 8 keV band after summing the MOS 1, MOS 2 and the 
PN data. Th e small number of cou nts in the nucleus in our Chandra 
observation (Blundell et al. 2006) precluded any detailed analysis. 
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Figure 1. Three colour plot of XMM-Newton data (188" X 116" image with 2 arcsec pixels; the image has been smoothed by a Gaussian with a of 2 pixels): 
0.2 — 1 keV band (red), 1 — 3keV band (green) and 3 — lOkeV band (blue) all on a square-root scale. The white overlaid contours are VLA A-array 
1.425- GHz radio contours (0.4, 1.7 and 7mjy/bcam), displaying the core and two lobes of the radio source. The blue source to the west of 6C 0905+39 is 
an obscured AGN (see section |3~3t . 



All spectral data described in the present paper were fitted over the 
0.5 - lOkeVband. 



3.1 The extended emission 



The extended emission is discussed in detail in lErlundetalJfeOOTh ; 
however, we briefly mention it here in the context of its effect as a 
potential contaminant of the core spectrum. Fig. [TJ clearly shows 
the extended emission from the bridge of 6C 0905+39 and that 
the number of X-ray photons increases towards the core. The ex- 
tended emission is not jet emission (since this would be very narrow 
and so would not be resolved by Chandra), but inverse-Compton 
up-scattering of the co smic microwave back ground by spent radio 
plasma from the lobes (Blundell et al. 2006). 

The bridge spectra (excluding a 15 arcsec region around the 
core) were extracted from the EPIC-pn, MOS 1 and MOS 2 data 
and were fitted with a Galactic-absorbed power-law giving a pho- 
ton index of V = 1.75 j^'^ an ^ a normalisation of 2.51q 3 x 
10~ 6 ct keV -1 cm -2 s _1 with a reduced-chisquared of xl = 0.9 
with 90 degrees of freedom (d.o.f.). The absorption-corrected X- 
ray luminosity of the extended emission in the 2 — 10 keV band is 



Lx = 1. 



j+0.3 



x 10 



44 erg s 1 



10 s 




partial co\ering model 
power law model 



3.2 The core 



Energy [keV] 

Figure 2. The spectrum of the core of 6C 0905+39. The solid line shows the 
best fit partial covering model and the dotted line shows the continuation of 
the directly observed power-law component. 



Only the PN data have been used in the following fits because each 
MOS spectrum contains roughly a quarter of the counts in the PN 
spectrum but a similar amount of noise. Spectral analysis shows 
that the MOS spectra are consistent with the PN spectrum. 

The XMM-Newton spectrum of the core is strikingly hard. 
A simple Galactic and intrinsic absorbed power-law is a poor 
fit with xl = 2.7 for 8 d.o.f. ( X 2 = 21.7). The best-fit pa- 
rameters give a photon index of T — 0.93, a normalisation of 
1.3 x 10 -6 ct keV -1 cm -2 s _1 and no intrinsic absorption. The 
spectrum (Fig. [2} appears to have a power-law component with 



an excess in the hard X-rays which can, in principle, be repro- 
duced by cold reflection (XSPEC model wabs*zwabs*[powerlaw 
+ pexrav]). This model is not a good fit however giving xl = 2.2 
for 7 d.o.f. (x 2 = 15.4) with an implausibly large amount of re- 
flection (reflection fraction ~ 33) and a photon index of T = 1.6. 

A partial-covering model with a redshifted iron line is a good 
fit CFig. [2j . Including an iron line, the model gives a covering frac- 



tion of 0.97±° ° 2 3 . 



intrinsic absorption of 2.6^J'° x 10 24 cm 



photon index T = 2.3^°' 3 and a normalisation of 5.0 



+0.7 



© 0000 RAS, MNRAS 000, 000-000 



3 



10~ 5 ct keV" 1 cm -2 s _1 for the power-law. The redshifted iron 
line (modelled as a Gaussian) was fixed at 6.4 keV in the rest- 
frame of the source with o = O.lkeV, giving an equivalent 
width EW ~ 0.9 keV; the normalisation of this component is 
2.2±\]l x 1CT 5 ct cm~ 2 s _1 . This gives xl = 0.9 with 6 d.o.f. 



5.4). Without an iron line, the fit values are similar: a cov- 



ering fraction of 0.95 



+0.01. 
0.04' 



an intrinsic absorption of 2.1 ._ 



10 cm ; a photon index Y — 2.1_ ' 2 and a normalisation of 



3.2 



+0.6 



x 10 _t> ct keV" 1 cm" 



" 1 , giving xl = 1-1 f° r 7 d.o.f. 



(X = 7.8). Including an iron line improves the fit by Ax = 2.4, 
so the data marginally favour its addition. The absorption-corrected 
luminosity in the 2 — 10 keV band for the partial covering model 



with an iron line is Lx 



3.2 



+0.3 



x 10 erg s and the ob- 



served flux is Fx = 3 x 10~ erg cm - s~ (the errors are 
poorly constrained). The column density of 6C 0905+39 is that 
of a Compton-thick AGN and its luminosity implies that it is a 
quasar. Such Comp ton-thick quasars are rare, Fig.[3](adapted from 
Gandhi et al. 2006) illustrates this point. It shows a selection of X- 
ray selected type-2 quasars which are not reflection dominated. We 
infer from this that 6C 0905+39 is particularly intrinsically lumi- 
nous and highly obscured. 

The paucity of counts in the Chandra data (14 ± 4 in the 
0.5 — 7keV band) means tha t no analysis was p erformed on the 
original detection of the core dBlundell et alT 2006). The Chandra 
data are consistent with the XMM-Newton data in as much as they 
are incapable of constraining the spectral shape of the nucleus. 
However, in the 2 arcsec region extracted for the Chandra analysis 
there are significantly fewer soft photons than in the XMM-Newton 
data (15 arcsec region). To quantify this we measure the observed 
count-rate below 1.5 keV, which is 1.5±0.9x 10~ 4 count s -1 , and 
compare it to that predicted from a model of the XMM-Newton data 
in the same band, which gives a count-rate of 5 x 10 -4 count s~ 1 . 
This deficit is significant at ~ 4a level and implies that the soft 
photons are from the extended emission which is less of a contam- 
inant in the smaller Chandra aperture. 

Given that there is a considerable amount of extended emis- 
sion from the lobes of the radio source (clearly seen in Fig. Q3 
that is likely to be contaminating the XMM-Newton nucleus spec- 
trum, we now try a highly obscured AGN with a power-law com- 
ponent at low energies from the extended emission. In order to test 
this, we fit a model consisting of Galactic absorption; a power- 
law with a photon index fixed to the value found for the ex- 
tended emission (r = 1.75, see Erl und et al. I l2007l) : an intrinsi- 
cally absorbed power-law (with a photon index fixed to T = 2 
as the data are unable to constrain it); and a redshifted iron line 
with it energy fixed at 6.4 keV and a — 0.1 (XSPEC model 
wabs[powerlaw + zwabs(powerlaw + zgauss)]). The best- 
fitting model is shown in Fig. [4] The intrinsic column density is 
3.5±o 4 x 10 2 ' 
2.7ln 4 x lO - 



cm and the normalisation of this component is 



' ct keV -1 cm' 



. The redshifted iron line has 



a normalisation of 2. 7^^^ x 10 4 ct cm 2 s 1 . This model gives 
a good fit to the data with xl = 1.1 for 7 d.o.f. (x 2 = 7.4); 
statistically the partial covering model is a slightly better fit with 
X 2 = 5.4 for 6 d.o.f. which gives A% 2 = 2. The observed 
flux is F x = 2.8to 4 x 10 _14 ergcm- 2 s _1 . The absorption- 
corrected luminosity in the 2 — 10 keV band of the central source is 
£*x,corc = 1.7to'i x 10 45 erg s _1 . The normalisation of the con- 
taminating power-law is LSj^'j x 10 -6 ct keV -1 cm -2 s _1 and 
its absorption-corrected X-ray luminosity in the 2 — 10 keV band 
is I/x.ext = l-llE) 2 x 10 44 ergs -1 . The luminosity of the ex- 
tended X-ray emission is Lx = 1.9l 'g x 10 44 erg s _1 implying 
that roughly a third of the power in the emission due to inverse- 
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Figure 3. This figure is adapted from lGandhi etal. (2006) (their Fig. 7, see 
their caption for details of how this plot was produced and the relevant ref- 
erences for each source). The X-ray luminosity axis shows the absorption- 
corrected luminosity in the 2 — 10 keV band of each source. B2 0902+39 
tB2_0902. lFabian et al.|[20o3) . SWIRE J104406.30+583954. 1 and SWIRE 
J1044 09.95+585224.8 (SW104406, SW104409 respectively; IPolletta et al] 
2006) and 6C 0905+39 have been added. This selection of X-ray selected 
type 2 quasars from the literature excludes sources which are clearly reflec- 
tion dominated. This implies that 6C 0905+39 is one of the most intrinsi- 
cally luminous and one of the most obscured quasars yet detected. 




Energy (keV) 

Figure 4. The best fit obscured AGN with contamination is shown in black. 
It consists of an absorbed Galactic and intrinsic absorbed power-law (green 
line); an additional power-law component fixed to have the same photon 
index as the extended emission (r = 1.75, red line) and a redshifted iron 
line (blue line). 



Compton scattering of the cosmic microwave background comes 
from close to the core. 
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T\ 



channel energy (keV) 

Figure 5. The PN spectrum of the second obscured AGN in the field (the 
blue source to the west of 6C 0905+39 in Fig.[T}. The solid line is the best 
fit absorbed power-law model with a redshifted iron line. 



3.3 Another obscured AGN in the field? 

The hard X-ray source to the west of 6C 0905+39 (the blue source 
in Fig. Q] located at RA 09h08m08.3s, Dec +39d43ml2s in J2000 
coordinates) is an SDSS source and so, following the SDSS nam- 
ing convention, is SDSS J090808.36+394313.6. It has a spectrum 
consistent with it being an obscured AGN (see Fig. [5} with column 
density of Nu ~ 1.28^^35 x 10 24 cm -2 when the power-law 
component is fixed at T — 2 because it is otherwise unconstrained 
(although varying it between reasonable values does not signifi- 
cantly alter the inferred value of the column density). Adding an 
iron line redshifted from 6.4 keV improves the fit from 1.4 for 4 
d.o.f. to 0.2 for 3 d.o.f., i.e. by A% 2 = 5. The redshift is mea- 



sured as z — 0.45, 



. There are a total of 261 counts (159 back- 



ground) in the 0.2—8 keV band, meaning that the derived quantities 
are somewhat uncertain; however, as Fig. Q] shows there are very 
few photons below 3 keV, this source must be a highly obscured 
AGN. The 2-10 keV band observed (absorbed) flux is F x = 
3. llti'o x 10 -14 erg cm -2 s" 1 and the absorption-corrected X-ray 
luminosity in the same band is Lx = 1.3^q'2 X 10 44 erg s _1 . 



4 DISCUSSION 

The spectral fitting of the nucleus of 6C 0905+39 presented in Sec- 
tion [3]2] shows that a model consisting of an obscured quasar nu- 
cleus with a contaminating power-law from the lobes is a good fit 
to the data. Its column density and luminosity clearly means that 
this is a Compton-th ick quasar. 

I Seymour et alj d2007t) have performed a Spitzer survey of 69 
high-redshift radio galaxies and find that the monochromatic lumi- 
nosity of 6C 0905+39 at 8 \im is Lg Mm ~ 5 x 10 44 erg s" 1 . They 
do not have mid-infrared data for 6C 0905+39, but in the majority 
of cases where they have longer wavelength data, the spectral en- 
ergy distributions (SED) of their high-redshift radio galaxies con- 
tinue to rise, often by one to two orders of magnitude. Observing 
6C 0905+39 at these wavelengths is important to confirm its nature 
as a Compton-thick quasar. 

Radio galaxies have long been supposed to be radio quasars 
oriented closer to the plane of the sky. Given its immense linear 
size of nearly 1 Mpc, 6C 0905+39 is likely to lie close to the plane 
of the sky and so, following the unification model, we are seeing 
the core through the torus and thus we would expect the source to 
be Compton-thick. This also means that beaming is unimportant. 



The measured obscuring column and absorption corrected lu- 
minosity suggests that 6C 0905+39 is a genuine Compton-thick 
quasar. (We note that the column density measured directly from 
X-ray spectr a depends on th e geometry of the obscuring medium 
as shown by iMatt et alJI 1999 . However, in the present case a sim- 
ple torus geometry seems relevant.) There are very few C ompton- 
thick quasars known at high redshift. iNorman et al. (2002) detected 
a good Compton-thick quasar candidate in the Chandra deep field 
South with 2 > 3. Apart from such clear examples, very few 
Compton-thick quasars are known especially at redshift z ~ 2, de- 
spite the fact that obscured AGN are expected to make up ~ 30 per 
cent of the hard X-ray background and 70 per cent should be found 
between 1 < z < 3, accor din g to populations synthesis models 
(see e.g. lNormanetalJl2002l and lGilli et al1l200ll) . 

The absorption-corrected luminosity in the 2 — 10 keV band 
ca n be converted to the bolome tric luminosity using results derived 
bv lVasudevan & Fabiarj 120071) who found that low and high Ed- 
dington ratio3 AGN have a bolometric correction of 19 ± 4 and 
55±13 respectively. In order to determine whether 6C 0905+39 is a 
low- or high-Eddington source, the black hole mas s is estimated us- 
ing th e black-hole-mass-bulge-mass relationship ( Magorria n et al.l 
1998) and the stellar mass (measured by ISevmour et al. 1 120071) as 
a proxy for the bulge mass. This gives a black hole mass for 
6C 0905+39 of Mesh ~ 2 x 10 9 M e . This is c onsistent with the 
value calculated using the method presented in Ma rconi & Huntl 
(2003), who make use of the black-hole- mass-bulge-luminosity 
relation, and the K-band magnitude (from ISevmour et al.l |2007|) 
along with a crude redshift correction. This assumes that such re- 
lationships hold to out to redshift z = 1.88 which is not nec- 
essarily true, but radio galaxies at lower redshi fts tend to be the 
hosts of similarly massive black holes (e.g. iLiu et al] l2006h so 
this estimate seems reasonable. Using these estimates for black 
hole mass, 6C 0905+39 appears to fall into the high-Eddington 
regime (I/boi/^Edd > 0.1) so the appropriate bolometric correc- 
tion is 55 ± 13 which implies that the bolometric luminosity of 
6C 0905+39 is L bol = 1.0 ± 0.4 x 10 47 erg s" 1 . This means that 
it is a very powerful source. Compared to other powerful, obscured 
quasars (see Fig. [3] assuming that the measured column density is 
accurate), 6C 0905+39 is amongst the most powerful, (X-ray de- 
fined) type 2 quasars currently known. 

Recent research has shown that high - Eddington objects tend 
not to be je tted (e.g. IPanessa et all 120071 : Isikora et af]|2007l and 
lMaozll2007h . However, 6C 0905+39 does not seem to be follow 
this trend as it is an extremely powerful source which is host to 
a radio source which spans 945 kpc on the sky and appears to still 
be actively fed. Its giant size indicates that it has been continuously 
active for at least the past ten million years. Such high-redshift giant 
radio galaxies may act as beacons for extremely powerful AGN. 



5 CONCLUSIONS 

Our XMM-Newton spectrum of 6C 0905+39, which is one of the 
highest redshift giant radio galaxies known, shows it to contain a 
Compton-thick quasar. 6C 0905+39 is one of very few such sources 
currently known and one of even fewer at high redshift. It is likely 
to be a high-Eddington source with a bolometric luminosity of ~ 
making it one of the most powerful sources known. 



10 47 erg s" 1 



1 where the Eddington ratio is defined as the bolometric luminosity divided 
by the Eddington luminosity 
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